Astronomy & Astrophysics manuscript no. aal8842-12 © ESO 2013 

April 29, 2013 



m 

o 

(N 

< 

(N 



o 






1 



3 



The progenitors of magnetic white dwarfs in open clusters 

B. Kiilebi^'^'^, J. Kalirai'*'^, S. Jordan^, and F. Euchner^ 

Institut de Ciencies de I'Espai (CSIC-IEEC), Facultat de Ciencies, Campus UAB, Torre C5-parell, 2" planta, 08193 Bellaterra, 

Spain 

Institute for Space Studies of Catalonia (ICE), c/Gran Capita 2-4, Edif. Nexus 104, 08034 Barcelona, Spain 

Astronomisches Rechen-Institut, Zentrum fUr Astronomie der Universitat Heidelberg, Monchhofstr 12-14, D-69120 Heidelberg, 

Germany 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore MD, 21042, USA 

Center for Astrophysical Sciences, Johns Hopkins University, Baltimore, MD, 21218, USA 

Swiss Seismological Service, ETH Zurich, Sonneggstrasse 5, 8092 Zurich, Switzerland 



Received XXXX, 2012; accepted XXXX, 2012 

ABSTRACT 

Context. White dwarfs are the final stages of stellar evolution for most stars in the galaxy and magnetic white dwarfs (MWDs) rep- 
resent at least ten percent of the whole sample. According to the fossil-field hypothesis magnetic fields are remnants of the previous 
stages of evolution. However, population synthesis calculations are unable to reproduce the MWD sample without binary interaction 
or inclusion of a population of progenitor with unobservable small-scale fields. 

Aims. One necessary ingredient in population synthesis is the initial-to-final-mass relation (IFMR) which describes the mass-loss 
processes during the stellar evolution. When white dwarfs are members of open clusters, their evolutionary histories can be assessed 
through the use of cluster properties. This enables an independent way of determining the mass of their progenitors. The discovery of 
the magnetic WD 0836+201 in the Praesepe cluster prompted the question whether magnetic fields affect the IFMR. In this work we 
investigate this suggestion through investigations of all three known MWDs in open clusters. 

Methods. We assess the cluster membership by correlating the proper-motion of MWDs with the cluster proper-motion and by an- 
alyzing the candidates spectroscopically with our magnetic model spectra in order to estimate the effective temperature and radii. 
Furthermore, we use mass-radius relations and evolutionary models to constrain the histories of the probable cluster members. 
Results. We identified SDSSJ085523. 87+ 164059.0 to be a proper-motion member of Praesepe. We also included the data of 
the formerly identified cluster members NGC 6819-8, WD 0836+201 and estimated the mass, cooling age and the progenitor 
masses of the three probable MWD members of open clusters. According to our analysis, the newly identified cluster member 
SDSSJ085523. 87+ 164059.0 is an ultra-massive MWD of mass 1.12 ± 0.11 M©. 

Conclusions. We increase the sample of MWDs with known progenitor masses to ten, with the rest of the data coming from the 
common proper motion binaries. Our investigations show that, when effects of the magnetic fields are included in the diagnostics, 
the estimated properties of these cluster MWDs do not show evidence for deviations from the IFMR. Furthermore we estimate the 
precision of the magnetic diagnostics which would be necessary to determine quantitatively whether magnetism has any effect on the 
mass-loss. 



cd 



Key words. Stars: white dwarfs - stars: magnetic fields - stars: distances - open clusters and associations: individual 

o 

C^ 1. Introduction to our understanding of the st ar formation efficiency in galaxies 

^_ ; dSomerville & Primack|[T999l) . 
^ ■ White dwarfs are the final products of stellar evolution for more 

• ^H ' than 95% of the stars in our galaxy. One of the most im- The fir s t clus t er IFMR studies we r e per formed by 

^ , portant tools for the understanding of stellar evolution is the IWeidemannI (Il977h : iKoester & WeidemannI (Il980l) . With now 

^ • initial-to-final-mass relation (IFMR), which, links the progen- about 40 white dwarfs associated with 11 clusters the empiri- 

itor m ass of a star to the final white dwarf mass d WeidemannI cal information on the IFMR has more than doubled o ver the 

Il977h . The IFM R is inherently connected to mass-loss mech- last few years (see [K alirai et al. 2009; Casewell et al. 120091 : 

anisms (see e.g. IWeidemann 2000). The understanding of stel- iDobbie e t al. 2009; IWilliams et al.i2009L for a recent discussion 

lar evolution from the perspective of the IFMR enables us to of cluster IFMR data). The increase in the data not only estab- 

constrain the critical mass necessary for Type II supernova ex- lished the support for the akeady estimated general trend, in 

plosions. Along with the initial mass function, this information which higher stellar masses yield more massive white dwarfs, 

can be used to estimate the birthrates and energetics of the su- but also resulted in an increase in the scatter. This spread is ar- 

pernovae, as well as the birthrates of the SN II products, the gued to be the result of differences in the host environments, 

neutron stars (Ivan den Bergh & TammannI Il99lb . The precise iMarigol (1200 ll) argued that one of these effects is metallicity, 

IFMR also constrains the galactic chemical evolution (i.e. en- for which metal-rich stars underg o heavier mass-loss and yield, 

richment in the interstellar medium), which in turn contributes on average, lighter white dwarfs. iKalirai et al.l (2005) was first 

to show this to be likely through spectrosco pic observations 
of white dwarfs in NGC 2099, and later Kal irai et all (|2007|) 
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Another physical process affecting the final mass of the 
white dwarf was proposed to be rotation. The fact that angu- 
lar momentum acts as an extra pressure support against gravity 
in side a degenerate structure has been known sinc e the works 
of lAnandl (Il965l) : iRoxburghl (IT965I) : lOsrike^ (IT966I) . Within the 
context of stellar evolution simple numerical methods show that 
this extra pressure causes the rotating degenerate core to expand 
and hence keep the m aximum temperature low er than the carbon 
ignition temperature (iDominguez et al.lll996h . This lower tem- 
perature enables the growth of the C-O core in mass while avoid- 
ing carbon ignition. In this way ultra-massive (1.2-1.4 Mq) 
white dwarfs can be formed from progenitors with masses be- 
tween 6-8 Mq. 

Magnetic fields, alongside angular momentum, were also 
proposed to be an effective factor of the mass-loss, after the dis- 
covery of a magnetic white dwaif (MWD) WD 08 36-1-201 as a 
probable member of the Pr aesepe cluster (EG59, Clav er et al.l 
|200U [Catalan et aP l2008al). Note that it was mislabeled in 
both works as explained in ICaseweU et al.l (l2009l) . Considering 
WD 0836-1-201 as a Praesepe member, and not taking the mag- 
netic field into consideration, the estimated white dwarf mass 
was higher than expected for its progenitor mass. If magnetism 
is an important factor affecting the IFMR, this would have im- 
plications for the MWD population. 

MWDs constitute at least 10% of the white dwarfs i f obser- 
vation al biases are considered (Liebert et al. 2()(B; Kawka et alJ 
l2007h . The current known population of MWDs has been in- 
creased considerably by the Sloan Digital Sk y Survey (SDSS) 
to about 220 objects p ansicke et al. 2002; Sch midt eta"l]l2003 : 
IVanlandingham et all BOOS; Kul ebi et alJ I2009). SDSS also dra- 
matical ly increased the tot al known white dwarf population (see 
e.g. Kle inman et alJ[2013L for Data Release 7) and recent stud- 
ies indicate that the number of M WDs in the SDSS could be 
as large as 521 (Kepler et al1 l2013h . The origin of the magnetic 
fields of MWDs is still under question. The accepted picture is 
the fossil-field hypothesis, where the field strengths are reminis- 
cent of an earlier stage of stellar evolution. Due to the conserva- 
tion of flux, the field strengths are amplified by the contraction of 
the stellar core to a white dwarf. For the case of MWDs, chemi- 
cally peculiar Ap a nd Bp stars were proposed as the progenitors 
(lAngeletal.lll981h . 

One problem of the fossil- field hypothesis is the relatively 
massive nature of the MWDs ('Lieberti ll988t IVenne s & Kawka 
l2008h . Whil e the mean v alue of the masses of the MWDs 
is 0.78 Mq (JKawka et alJ 12007). the mean m ass of the non- 
magn etic white dwarf sample is 0.661 Mq (iGianninas et alJ 
1201 lb . The hypothesis that this could be a result due to the 
impact of m agnetism on the ma ss-loss was considered by 
IWickram asinghe & Ferrariol (l2005l) via population synthesis. 
Their conclusion was that the current number distribution and 
masses of high-field magnetic white dwarfs (HFMWDs, B > 
10^ G) are not generated by an inclusion of a modified IFMR, 
but rather by considering the contribution of about 10 percent of 
A/B stars with unobservable small scale magnetic fields. 

The aim of this work is to analyze the MWDs in clusters with 
sophisticated magnetic models, in order to estimate the progen- 
itor masses. This is undertaken by calculating the masses of the 
MWDs and their progenitors by using the additional information 
estimated due to their membership in open clusters, namely their 
distances and total evolutionary ages 

Up to now, only two MWDs have been discovered as mem- 
bers of open clusters, one in Praesepe ( WD0836-H201) an d 
the other one in NGC 6819 (NGC 6819-8. lKaiirai et al.ll2008h . 
Neither of them has been analyzed with magnetic models so 



far In order to increase the statistical sample further, we cor- 
relate the properties of SDSS white dwarfs that show evidence 
of magnetic fields with the kinematic properties of open clus- 
ters in order to determine the possibility of their cluster mem- 
bership. While we identified only SDSS J085523. 87+ 164059.0 
(SDSS J0855-I-1640) as a member of an open cluster, we also 
analyzed NGC 68 19-8. Furthermore, in this work we discuss for- 
mer investigations on the influence of magnetism on the IFMR, 
compare our results with previous studies, and evaluate the anal- 
ysis employed. 

In Sect. |2] we assess the cluster membership of the currently 
known SDSS MWDs in the literature through proper-motion in- 
formation. Afterwards, we describe our methods to analyze the 
cluster MWDs in Sect.[3]and explain what the relevant uncertain- 
ties are. In Sect.|4]we apply our method to cluster MWDs and in 
Sect.|5]we discuss whether we can conclude on possible eff'ects 
of the magnetism on the IFMR 

2. Cluster Membership 

In our investigation we compared the astro metric properties of 
137 SDSS hydrogen-rich (DA) MWDs (see iKiilebi et alj|2009l 
and references therein) with the properties of 520 open clus- 
ters fr om the Catalogue of Open Cluster Data (iKharchenko et al.l 
|2005|) . In order to assess the membership probabilities of the 
white dwarfs based on their k inematics, we use the proper- 
motions estimated by the SDSS dMunn et alj|2004 . The proper- 
motions of the MWDs were acquired from the CAS JOBS SQL 
interface (Li & Thakar 2008). 

Our criterion for open cluster membership was through 
checking whether the MWD lie within the tidal radius of an open 
cluster (values taken from the catalogue of lPiskunov et aLEOOSh 
and whether the proper-motion of the MWD is within 3 cr of 
the cluster proper-motion. As a result we obtained only one new 
possible open cluster member SDSS J0855 + 1640 of the open 
cluster NGC 2632 (Praesepe), which is a nearby open cluster 
that has b een the subject of many white dwarf and IFMR inves- 
tigations (Luvten 1962:'Egg en & Greensteinl[l965l:IClaver et al.l 
2001). Additionally with our calculations we recover the proto- 
typical cluster member MWD, WD 0836-1-201. 

Another constraint on possible membership is the position 
of the white dwarfs in the color-magnitude diagram (CMD) with 
resp ect to th e cooling sequence expected for the cluster (see 
Kalir ai et alJ 12008). Since this method involves a prescription 
for the parameters we want to evaluate, e.g. deviation from the 
expected IFMR due to the effect of magnetism, we did not un- 
dertake such an analysis. Rather, the evolutionary constraints on 
the cluster membership come from the estimated cooling age of 
the white dwarf, i.e. the cooling age cannot be larger than the 
cluster age. 

Finally the membership possibility of both 
SDSS J0855-H1640 and WD 0836-H201 can be assessed through 
considering the studies made on the WD populations within 
the open clusters. Both of these MWDs are possible members 
of Prasepe. The current number of white dwarfs in this cluster, 
confirmed by radial velocity observations of 'Casewe ll et al.l 
(2009) is eight, including the known magnetic WD 0836-1-201. 
compare d to numerical simula t ions o f the dynamics of open 
clusters dPortegies Zwart et al.l 12001 1) the number of white 
dwarfs found in Praesepe is lower than expected, especially if 
compared to the number of observed giants. Hence it is logical 
to expect more white dwarfs to be discovered if the search is 
extended at least up to the tidal radius, as has been done in our 
work. 
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Table 1. Photometric and astrometric properties of the DAHs with possible Praesepe memberships. 



MWD 



RA / (h m s) Dec. / (deg ' ") r / mag 



'■ / mag 



: / mag RA p.m. / (mas/yr) Dec. p.m. / (mas/yr) 



WD 0836+201 
SDSSJ0855+1640 



08 39 45.56 
08 55 23.87 



+20 00 15.7 
+ 16 40 59.0 



18.11= 
18.80= 



0.01 
0.01 



18.36±0.01 
19.05±0.02 



18.66= 
19.32= 



;0.04 

0.08 



-32.93±2.87 
-33.14±2.75 



-15.98H 
-14.71H 



:2.87 
:2.75 



Notes. The columns indicate the right ascension (RA), declination (Dec); the SDSS photometric magnitudes r, i, z which are relevant for this 
work; and finally the proper-motions. Compare the proper-motions of the objects with the proper-motions of Praesepe RA p.m.= -35.90 ± 
0.13 mas/yr, Dec. p.m. = -12.88 ± 0.1 1 mas/yr ( iKharchenko et al. 200^ . 



3. Analysis 

One of the necessary ingredients to determine the cooling age 
of a white dwaif is the determination of T^^ and \ogg values. 
Although this is undertaken by x^ fits to the spectral lines for the 
non-magnetic white dwarfs, the situation for the MWDs is more 
complicated due to the influence of magnetism on the atomic 
transitions and the distribution of the magnetic field strengths 
over their surface. 

To overcome this difficulty and to analyze the magnetic field 
geometry of the MWDs, we use a two-step analysis approach 
in which first the model spectra for magnetized atmospheres are 
calculated for a distribution of magnetic field vectors with re- 
spe ct to the fine of sight and the n ormal on the surface of the star 
(see lJordanlll992t|jordan & Schmidt.2003.) . This is repeated for 
given Teff and \ogg values and a database is generated. 

Later these spectra representing individual vectors could be 
accessed from the database and added up in order to calculate 
the spectrum for a given geometry. The parameters of the field 
geo metry are determine d by an evolutionary fitting procedure 
(see lEuchner et al.ll2002l) . 

This approach was used by iKtilebi et alJ (l2009h where the 
models in terms of offset dipoles for the single phase of these 
white dwarfs were determined. However, for that work the T^^ 
values were determined by the photometric colors. In this current 
work, we use the same database and procedure to determine the 
Teff value of the MWDs in detail. 

Unfortunately, for MWDs it is not possible to determine the 
log g values due to the difficulties in accounting for the simulta- 
neous impact of Stark and Zeeman effects on the spectral lines. 
Moreover, no atomic data exists for hydrogen in the presence of 
both magnetic and electric fields for arbitrary strengths and ar- 
bitrary angles between the two fields. Up to now, only simple 
cases of parallel e lectric versus magne tic configurations have 
been investigated dFriedrich et aH 1 19941) . Since in our models 
only crude approximations have been used, it is not possible 
to avoid the systematic uncertainties, especially in the low-field 
regime (< 5 MG) where the influence of the Stark effect on the 
spectral lines is strong relative to the effect of Zeeman split- 
ting. Therefore, we do not determine log g values from our spec- 
tral fits, but rather assume \ogg value of 8.0 for our magnetic 
model spectra database. Thi s approach have been apphed mul- 
tiple times in the literature (iGansicke et al.l 120021; iGirven et alj 
12010 ). Additionally, we estimated the error budget introduced 
by the undetermined surface gravity through repeating our fitting 
procedure for the cases in which log g -1 .Q and log g - 9.0, and 
later propagating the resulting uncertainties. 

The determined parameters are listed in Table|2] the errors 
in Teff are defined by the spacing of our grid points in our 
database. The observed spectra are compared to theoretical mod- 
els in Fig.[I] The fainter flux of SDSS J0855+1640 with respect 
to the flux of WD 0836+201 despite its higher effective temper- 
ature hints a smaller radius, hence a more massive nature, if we 
assume both MWDs to be at the same distance. 




5000 5500 6000 6500 



Fig.l. Spectral fits to WD 0836+201 (top) and 
SDSSJ0855+1640 (bottom) determined in this work. The 
parameters are given in Table|2] The inconsistencies in the 
WD 0836+201 line fits are due to the simple treatment of the 
relatively large Stark effect in the low-field regime (B < 5 MG). 



In order to determine the MWD masses we used the syn- 
thetic magnitudes (in the ugriz photometric system of SDSS) 
for carbon-oxygen (CO) core white-dwarf cooling models with 
thick hydroge n layer (M n/Mt - lO"'*) masses ( Fontaine et al] 
2001; Holber g & Berger on 2006 n; since these are non-magnetic 
models the effect of the magnetism is not included in these cal- 
culations. 

r, / and z-bands sample mostly the continuum hence these are 
the bands least affected by the Zeeman splitting of the Balmer 
lines (see Fig.|2]i. Moreover, in order to minimise the uncertain- 
ties of our analysis we restricted ourselves to the use of r, i and 
z. If we compare the estimated MWD mass from the / band with 
the masses derived from the u band we obtain a difference as 
high as 8%; in the case of masses estimated from / and g mag- 
nitudes, the differences are as small as 1%. The results for z and 
r magnitudes are consistent within 0.05%. This shows that the 
photometric magnitudes that correspond to spectral continuum 
are consistent with each other, thus we use the / band in our cal- 
culations. 

Given the photometric magnitudes and the estimated Teft 
of an MWD, we interpolated the synthetic magnitude grids of 
Holberg & Bergeron (2006) in order to estimate its mass and 
cooling age. The systematic eiTors depend on observation errors 
of the photometry and the eiTors of Teff from the fits. 

Furthermore, we subtracted the estimated cooling ages from 
the total age of the open cluster to calculate the progenitor ages 
of the MWDs. With this information at hand, estimating the ini- 



http : //www . astro . umontreal . ca/~bergeron/CoolingModels 



B. Kiilebi et al.: The progenitors of MWDs in open clusters 
Table 2. Model fits to the SDSS observations of the objects considered in this work. 



MWD 


Plate-MJD-FiberlD 


Sdipoic / MG 


Zoff / ''WD 


(■ / deg 


Tefl- / K 


WD 0836+201 
SDSSJ0855+1640 


2277-53705-484 
2431-53818-522 


2.83+0.19(0.62) 
12.6+1.0(3.9) 


-0.26+0.02(0.11) 
-0.25+0.02(0.16) 


26+4 (33) 
44+6 (40) 


17000+500 
20000+500 


NGC6819-8 


- 


10.3+1.1 (2.4) 


-0.20+0.03 (0.09) 


52+9 (22) 


19000+1000 



Notes. The first two columns indicate the SDSS name of the object; the plate, Modified Julian Date and fiber ids of the observations; the 
remaining columns indicate the model fit parameters for offset dipole models. The model parameter ; refers to the inclination of the magnetic 
dipole axis with respect to the line of sight, and the offset is along the axis of the magnetic field in terms of the stellar radius. The errors in 
parentheses refer to the case in which log g is not fixed to the value of 8.0, but is allowed to vary between between the values of 7.0 and 9.0. 




10000 11000 



Fig. 2. The diagram above shows the response curves of the 
ugriz photometric pass-bands with 1.3 airmass correction ap- 
plied, taken from the SDSS websit^l compared to the MWD 
spectrum of two different effective temperatures and same mag- 
netic field geometry which for this case is the 12.61 MG offset 
dipole model of SDSS J0855+1640. The fluxes are arbitrarily 
scaled, but consistent relative to each other. The effective tem- 
peratures are 20 000 K (top) and 15 000 K (bottom). The u, g 
and r bands cover the wavelengths which are effected by the 
magnetism in the MWD spectra for field strengths smaller than 
50 MG. However this effect is minor for the r band. 



tial masses is rather straightforward given the stellar evo lution 
models. Hence, we interpolated the iBertelli et al.l (12009 ) stel- 
lar evolution grids with cubic splines to calculate the progenitor 
masses given the progenitor lifetime, until the beginning of the 
Asymptotic Giant Branch phase. 



4. Results 

4.1. WD 0836+201 

Previous discussions of the pr operties of WD083 6-H201 are 
mostly based on the analysis by IClaver et al.l (1200 ih where the 
n components of the Zeeman splitted lines were used to deter- 
mine Fe ff ^ 17 098 + 350 K and log^ = 8.32 + 0.05. Using their 
values. .Catalan et aP (l2008a) re-estimated the mass and the pro- 
genitor mass o f this object to discuss the effect of the magnetic 
field on IFMR. ICasewell et alj <2009l) mentioned that the analy- 
sis with non-magnetic model spectra may lead to uncertainties in 
the temperature determination. From the SDSS colors (g-i, u-r) 
of WD 0836-1-201, they determined the effective temperature to 
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be about 15 000 K, and argued on this basis that the Teg and \ogg 
determination from the spectroscopic non-magnetic model fits 
are likely to be less robust. With our theoretical spectra for mag- 
netised white dwarf atmospheres we obtain a more reliable fit of 
Tsff = 17 000 + 500 K within the spectral modeling uncertain- 
ties, a result which is consistent with Claveret al. (2001). Since 
the estimation of the MWD mass relies on the cluster distance 
of Praesepe, additional errors come from the distance errors. 

There are multiple distance measurements of Praesepe. Here 
we adopt the value of 184.5 + 6 pc as the d istance to the 
cluster center, which was used by ICasewell et al. (2009). It is 
the weighted mean of the Hippa rcos-based distan c e mea sure- 
ment (m - M)o = 6.24 + 0.12 iMermilhod et all (I1997D . the 
ground-based parallax meas urement (m - M)o = 6.42 + 0.33 
(iGatewood & d e Jonge'1994') and the photometric determination 
(m - M)o = 6.33 + 0.04 (An et al. 2007). It is not only a consis- 
tent approach to account for the multiple distance estimations in 
the literature, but also advantageous for us since we would lik e 
to compare our estimates to the results of lCasewell et al.l (l2009l) . 

For the determination of the progen itor ages, we rehed 
on the cluster age of 625 + 50 Myr (.Perrvman et al.l 1 19981 
from isochrone fits with a metallicity of Z = 0.024). We 
have chosen a value of Z = 0.027 taken from ICatalan et al.l 
(2008a) but our mass and age determinations do only weakly 
depend on the choice of metallicity: If we instead use the low- 
est recent value for the metallicity of Praesepe (Z = 0.019, 
IClaveretal.]l200U ICaseweU et"ani2009l) flie progenitor age be- 
comes only 1% longer, which implies the progenitor mass esti- 
mate s would be 0.03 Mfr^ s maller, consistent with the estimates 
from lCasewell et all (l2009l) . 

Our results are summarized in Table|3] It should be noted that 
the smaller errors of the former works on WD 0836 -f 201 are due 
to underestimation of magnetic effects in spectral fits. Moreover, 
we can only constrain the status of WD 0836-1-201 as an outlier 
to the semi-empirical IFMR up to a limited confidence interval. 
Quantifying this depends strongly on the uncertainties in the Teif 
errors, which are adopted as the difference between the best fit 
grids in our models. However, when a greater level of confidence 
is applied to the observables used for the parameter estimation, 
the resulting values for WD 0836-1-201 entail the semi-empirical 
IFMR, indicating no evidence for the possible effect of mag- 
netic fields on the mass-loss. This situation is depicted in Fig.|4] 
where we use 3 cr errors for the photometric magnitudes and use 
2 000 K asa more realistic error for the Tetf . 

4.2. SDSS J0855+1 640 

There are certain caveats for the membership status of 
SDSS J0855-I-1640 when the rest of the Praesepe members are 
considered, since it is further away from the cluster center and 
lies at the limb of the circular region defined by the tidal radius 
of the cluster Within this radius the cluster members are ex- 
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Table 3. log g, masses, radii and cooling ages estimated with the photometric analysis. 



MWD 



log g I dex Mwd/ Mq radius / 0.01 Rq /cool / Myr 



,/Myr MJ Mq 



WD0836+201 8.33 + 0.08 0.82 + 0.05 1.02 + 0.06 

SDSSJ0855+1640 8.84 + 0.22 1.12 + 0.11 0.67+0.14 



221+" 



391+ 



my% 



234+177 



Q 144+0.181 



1 Sr 8+1-870 



NGC6819-8 



7.77 + 0.10 0.50 + 0.05 



1.53+0.76 



66!!!! 



243011° 



1 570+"°^'^'' 



pected to have the same kinematic properties as the cluster and 
high cluster proper- motions relative to the field stars enable us 
to discriminate the members. 

The value for the tidal radius of Praesepe differs in the lit- 
erature, and is incre asing as the studies be come more recent. 
Starting with 12 pc (Mermilliod et al. 1990'), a higher value of 
16 pc was obtained bv Adams et al.. (.2001.) . an d later an even 
higher value of 17.1 ± 1.2 pc was estimated bv iPiskunov et alJ 
( l2007h . In our wor k we used the latest value of 18.6 + 1.4 pc 
value of lPiskunov 'et al. (2008). 

If we assume that SDSS J0855+1640 lies at 184.5 pc, which 
is the distance to the core of the cluster, its angular separation 
from the cluster center of 4.58° translates into a spatial distance 
of 14.7 pc. Although numerous values exist in the literature, the 
distance of SDSS J0855-I- 1640 from the core puts it safely within 
the tidal radi us of Praesepe for all det erminations except the one 
estimated bv lMermilliod et al.l (Il990l) . 

Another support for the Praesepe membership of 
SDSSJ0855+1640 is the determined cooling ages. The 
value estimated through the T^s value of 20 000 K (see Table[3]) 
yields cooling ages smaller than the age of Praesepe. Hence the 
cluster membership of the SDSS J0855H-1640 cannot be out- 
ruled by evolutionary arguments. However, it is still necessary 
to carefully consider the spatial position and the distance of 
SDSS J0855+1640 from the core of the cluster 

Due to the mass segregation within a cluster as time pro- 
gresses, relatively massive stars (> 1 Mq) concentrate within 
the core. Since these stars evolve to white dwarfs, their popula- 
tion is also expected to be found within the core. This is con- 
firmed for the known white dwarfs of Praesepe and other clus- 
ters. However, there is a selection effect since spectroscopic in- 
vestigations specifically aim at the core of the cluster, due to the 
above reasoning. 

However, it is known that the number of observed white 
dwarfs inside open cl usters is lower than expected (Weidemann 
1 1 977t IWiniamsll2002h . The fact that white dwaifs tend to evap- 
orate from the inner core of an open cluster was already dis- 
cussed for the case of the missing white dwarfs in Hyades 
(IWeidemann et al.l 11992*). In this scenario, white dwarfs may 
leave the cluster by gaining small velocities (0.1 km/s) with re- 
spect to the cluster center Recent observations show that there 
is evide nce for this suggestion for the case of the Hyades open 
cluster (ISchilbach & Roeseiil201 Ih . The authors reason that, af- 
ter the formation of white dwarfs near the cluster center, due to 
their lower masses they are subjected also to mass segregation 
and move outwards. 

Another possibility for gaining a kick velocity of sev- 
eral k m/s is due to asymm etric mass-loss during the AGB 
phase (iFellhauer et al.l 120031) . This was tested observationally 
for the globular cluster NGC 6397 and found to be plausible 



(iDavis et al.ll2"008h . Since SDSS J0855+1640 is within the tidal 
radius we did not consider this case for the following assessment. 

Given the cooUng age of SDSS J0855+1640 of 390 Myr, a 
kick velocity of 0.1 km/s implies the distance travelled since its 
formation would be about 40 pc. This value is significantly larger 
than the tidal radius of Praesepe, and may explain the distance of 
SDSS J0855-I-1640 with respect to the cluster center Although a 
kick velocity explains the distance of SDSS J0855+1640 from 
the cluster core, at the same time does not guarantee its radial 
distance to be same as the cluster distance. If the MWD has 
gained a kick velocity in the radial direction, it may lie consid- 
erably behind or in front of the cluster center 

If the maximum distance travelled is approximately 40 par- 
sees, and the distance of SDSS J0855+1640 to the center of the 
cluster is 14.7 parsecs, then the maximum possible unaccounted 
radial distance equals to 37.2 parsecs. This is two times the tidal 
radius if the value of 18.6 parsecs is assumed for the value of the 
tidal radius. 

In order to estimate the maximum possible uncertainty, we 
calculated the progenitor properties of SDSS J0855 -1-1640 using 
two tidal radii as the distance error Within these uncertainties 
the inferred cooling age of SDSS J0855+1640 is smaller than the 
age of the cluster, hence its membership could not be outruled. 

Our analysis shows that SDSS J0855-(-1640 is possibly a part 
of the rare group of ultra-massive white dwarfs which have 
masses higher than 1.1 Mq. This group of objects are proposed 
to be final products of binary mergers and were observed nu- 
merously in EUV surveys (Bergeron et al. 1991; Marsh et al] 
1 19971) . Given the massive nature and the position on the IFMR 
of SDSS J0855+1640, a merger scenario could be probable. 
However, when a larger confidence interval is applied for cal- 
culating the errors of the initial and final mass of this star, these 
values would entail the semi-empirical IFMR (see Fig.lUi and a 
single star origin cannot be outruled. The nature of this star can 
only be understood conclusively, if its Teft errors could be con- 
strained to a value smaller than 2 000 K within 3 a confidence 
intervals. 

4.3. NGC 6819-8 



The spectrum of NGC 6819-8 is taken from lKalirai et al.l (l2008l) . 
where it was not analyzed in detail due to the faintness of the 
star and the low signal-to-noise. We have assumed the object to 
be a member, whereas in Kalirai et al. (2008) the membership 
assessment was made by comparing the observed V magnitudes 
with the theoretical magnitude from fitting the Balmer lines. The 
comparison indicated that these magnitudes were inconsistent by 
more than 1 cr, so the white dwarf was assumed not to be a part 
of the cluster Kalirai et al. (2008) used a strict cut to ensure that 
only real members are considered, however we cannot employ 
such a method for our MWDs due to the complications involved 
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Fig. 3. Spectral fit to NGC 6819-8. Fit parameters are given in 
Tabled 



in the modeling of the photometric magnitudes under the influ- 
ence of the magnetic fields (see Sect.O. We assumed the mem- 
bership of this object, based solely on the fact that its location 
agrees with the WD cooling sequence of the cluster 

Our spectral analysis yielded the Teft and the mass of the 
white dwarf. Subsequently we used the same methods used for 
the Praesepe white dwarfs to estimate the initial mass of the ob- 
ject. One difference of our analysis was our usage of V mag- 
nitudes rather than the / band. This might introduce extra errors 
since the V band is influenced by the splitting of hydrogen/? line. 
However, the estimated mass error should be smaller than 1%, 
enumerated in the comparison between estimates with / band 
which represents the continuum flux and the g band which cov- 
ers multiple Balmer lines (see Fig.|2]l. For the fits an additional 
error is introduced due to the short wavelength interval of the 
spectrum, which makes it harder to distinguish between models 
with different Teff . 

Due to the aforementioned reasons, the main uncertainty in 
the (final) mass of the NGC 6819-8 comes from the inaccuracy 
of the Tetf determination. However the error in the initial mass 
of the NGC 68 19-8 is not affected by the uncertainties of the 
spectral modeling since the error of the cluster age dominates 
over other uncertainties in the spectral analysis as well as the 
white dwarf cooling age itself. This means that the error of the 
progenitor mass of the object is comparable to the errors of the 
progenitor masses of other objects in the cluster 

The estimation of the cluster age is achieved by theoreti- 
cal isochrone fitting, which is dependent precisely on the dis- 
tance, metallicity and reddenning. The distance and metallicity 
for the cluster NGC 68 19 have been estimated through main se- 
quence fitting by Kalirai et al.l (1200 ll) . In their work the adopted 
reddening of E(B-V) = 0.10-0.14 yielded slightly larger values 
((m - M),. = 12.30 + 0.12) with respect to the previous studies. 
This was due to the reddening values being smaller Using these 
metallicity, reddening and distance p arameters, the isochrone s 
yield 2.5 Gyr for the age of the cluster (iKahrai et al.l200ill2008h . 
It is noted that the degeneracy of the input parameters might in- 
duce a about 10 percent change in the result — which is the error 
adopted for the age in this work — however, the parameters are 
in agreement with the literature and, more importantly the lower 
main sequence were reproduced well by their model. 
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Fig. 4. IFMR and the posit ion of the Praesepe white dwarfs (cir- 
cles, Casewefletal. 2009) and NGC 6819 white dwarfs (dia- 
monds, Kalira i et al.ii 2008) are shown in this diagram. The solid 
line is the best fit oflCa sewell et al. (2009), and the dashed line 
is the best fit from Kalirai et al. ( 2008) which is sensitive to the 
low-mass end. Red circles denote the positions of the MWDs 
discussed in this paper, and WD 0836+197 is a known outlier 
For the MWDs considered in this work, the bold (black) errors 
correspond to the estimations made with the standard procedure 
as explained in Sect.[3j whereas the thin (red) error bars corre- 
spond to the 3 (J error estimations including an adopted T^s error 
of 2 000 K. 



Using the V magnitude of 23.3 for NGC 6819-8, the distance 
modulus for NGC 6819 as 1 2.30+0.12, and a me taUicity of Z = 
0.017 (all parameters from iKalirai et al.l 120081) . we reached a 
progenitor mass of 1 .57 Mq which is similar to the rest of the 
white dwarfs in the cluster For comparison, we have applied 
our method to get the progenitor mass of the non-ma gnetic white 
dwarfs in NGC 6819, and reproduced the results of .Kahrai et al.l 
(l2008l) . 

Given the uncertainties in the cluster age, the estimated pro- 
genitor mass of NGC 6819-8 is the most precise for the known 
MWD population. Although the progenitor mass is relatively 
low, it is still within the mass range of Ap stars, which is the 
progenitor population suggested by the fossil-field hypothesis. 

After evaluating the errors and estimating the values of fi- 
nal and progenitor mass of NGC 6819-8, we end up with values 
comparable to the rest of the NGC 6819 white dwarfs. Within 
the uncertainties we conclude that there is no evidence of any 
effect of magnetism on the mass-loss history of NGC 6819-8. 



5. Discussion and Conclusions 

In this work we investigated the evolutionary histories of MWDs 
which are probable members of open clusters. First of all, we 
correlated the proper-motions of the known magnetic DA white 
dwarfs in the SDSS with the open clusters. Afterwards we esti- 
mated the properties of probable cluster MWDs with our mag- 
netic spectral modeling, and used the cluster properties to fur- 
ther constrain their evolutionary history. This method was ap- 
pHed to two currently known cluster MWDs (WD 0836+201 and 
NGC 68 19-8) for the first time and additionally to a newly dis- 
covered SDSS J0855+1640. 
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We also compared our results with the the former analy- 
ses which have been undertaken with the use of non-magnetic 
WD spectral fits. By that, we assessed the conditions of appli- 
cability of non-magnetic analysis for MWDs. Our results show 
that SDSS J0855+1640 is probably an ultra-massive MWD and 
NGC 6819-8 has properties very similar to the rest of NGC 6819 
members. 

Furthermore, through the use of cluster properties we also 
estimated the progenitor masses of these MWDs. In Fig.|4] we 
compare our results and the data from Praesepe and NGC 6819 
cluster white dwarfs within the scope of IFMR. Our calculations 
show that within the uncertainties of the magnetic analysis, all 
of the MWDs follow the general trend of the IFMR for non- 
magnetic white dwarfs derived from seven open clusters and 
measurements from 41 white dwarfs (ICasewell et al.l 120091) . It 
should also be noted that this general trend is not followed by the 
known outlier WD 0836+197, for which neither magnetic field 
nor rapid rotation was detected; it was proposed to be a pos- 
sible radial velocity v ariable or the product of a blue straggler 
(ICasewell et alJl2009h . 

The most important diagnostic tool in the estimation of the 
white dwarf and the progenitor masses is the determination of 
the effective temperature. The current progress in the magnetic 
spectral analysis allow for the estimating the best fit for this pa- 
rameter however the formal determination of the confidence in- 
tervals is still not possible due to the complications caused by 
the effect of magnetism on the spectra. Hence in our work we 
adopt the distance between the bins in our database as errors. 
Furthermore, when we repeat our calculations with 3 cr errors 
for Teff, we show that the the estimated progenitor and white 
dwarf masses would entail the semi-empirical IFMR values (see 
Fig.|4|. According to these calculations, to be able to conclu- 
sively determine the effect of magnetism on the mass-loss, the 
magnetic diagnostics should be able to attain Teff errors smaller 
than 2 000 K within 3 cr confidence intervals. In order to conclu- 
sively quantify the confidence levels of complicated ;t'^ topogra- 
phies, a Markov Chain Monte Carlo type of analysis is needed 
however, this is beyond the scope of this work. 

For the the high-mass white dwarf SDSS J0855+1640 we re- 
solve that it could also be the re sult of binary evolution as dis- 
cussed bv lBerperon et al.l (Il991l) and lMarsh et al.l (1 19971) . In this 
case it cannot be compared to the IFMR of single stars. The inci- 
dence of magnetism among ultra-massive white dwarfs is rather 
high (jx^nnes & Kawka 2008), suggesting a possible relation of 
the evolution histories, masses and the magnetic fields. 

The membership of SDSS J0855+ 1640 could be clarified 
through the de termination of its trigonometric p arallax e.g., with 
Gaia (see e.g. iTorres et al.l 120051: lJordanll2007l) . This measure- 
ment is also important for determining the origin of this star, 
since our mass determination is based on the assumption that the 
star is at the distance of the center of the Praesepe open cluster. 

For the low mass end, NGC 6819-8 presents similar proper- 
ties - if not a slightly lower progenitor mass - with respect to 
the rest of the NGC 68 19 cluster members. When the position of 
NGC 6819-8 in the IFMR plot is compared to the solid fine in 
Fig.JH which is the best fit of the IFMR for the low-mass end 
JKalirai et al. 2008), no considerable difference can be observed 
within the error bars. 

One should also note that there is considerable scatter for 
the estimated initial masses smaller th an about 3 M(?) wh ich was 
also observed for NGC 2099 (M 37, Kali rai et al.ll2005h . It was 
proposed that the theoretical uncertainties of mass-loss in this 
regime might also be the reason for this discrepancy. At this 
stage it is difficult to distinguish between the observed scatter in 



the low mass end due to other physical mechanisms (i.e. metal- 
licity), and the effect of magnetism. 

It is also possible to test the IFMR using double white 
dwarfs in wide binary systems, since these systems are 
assumed to have evolved without interacting. The num- 
ber of known MWDs in WD-WD wide b i naries is six, 
namely REJ 0317-853 (Barstow et al." '1995"; 'Ferrari o et aH 
1997), SDSS J 1300+5904 (Girvenetal . 2010) and re- 
cently SDSS J092646.88+ 132 134.5 (SDSS J0926+ 1321), 
SD SS J15074 6. 80+52 0958.0 (SDSSJ1507+5209, 

Dob bie et alJ lIoTl . SDSS J074853.07+ 3 02543 .5 and 
SDSSJ150813.24+394504.9 (Dobbie etdj I20TI . Their 
evolutionary status presents different challenges. 

The ultra-massive and hot nature of RE J 0317-853 with re- 
spect to its non-magnetic counterpart turns out to be problem- 
atic for a single star. This problem of "age dilemma" was re- 
solved by invoking a scenario where the system initially con- 
sists of three compo nents and REJ 0317 -853 is the result of a 
close binary merger (iFerrario et al.lll997l l. A similar evolution- 
ary dilemma was also observed for SDSS J1507 +5209, where 
the MWD component of the binar y is relatively hot a nd young 
compared to the non-magnetic one (iDobbie et al.ll2012l) . hence a 
triple system was also suggested as an explanation. The remain- 
ing MWD in wide WD-WD pairs do not show similar problems. 

A similar type of analysis can be undertaken with com- 
mon proper motion (CPM) binaries consisting a white dwarf 
and a FGK type star, where the total age and the metallicity 
of t he white dwa r f can b e assessed by analyzing its compan- 
ion. |Catallnetal] (12008B) undertook the first of these studies in 
which one of the objects (40EriB, WD 0413-017) is an MWD. 
Their observational data show dispersion with respect to the 
semi-empirical IFMR, however whether the cause of this dis- 
persion is magnetism and/or rotation has not been discussed in 
detail. It should be noted that the magnetic field strength of 40 
Eri B (^ 2.3 kG, Fabrika et al. 2000) is significantly lower than 
the magnetic field of MWDs in this work (> 1 MG). 

CuiTently there are ten MWDs with determined white 
dwarf and progenitor masses. Not considering the two objects 
REJ 0317-853 and SDSS J1507 +5209 which are believed to 
have formed by mergers, it has been suggested that magnetism 
could be a factor in m odifying the mass-l oss based on the discov- 
ery of WD 0836+201 (IClaveretal.l2001h . However our analysis, 
which also entails WD 0836+201, shows that when the effect of 
magnetism on the observed spectrum is included in the uncer- 
tainties, the deviation of the properties of MWDs from the semi- 
empirical IFMR disappears. The conclusive test for the extent of 
the dispersion caused in the IFMR by the magnetic fields will be 
possible through better magnetic diagnostics and more precise 
distance measurements. 
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